Holographic Charge Density Waves 
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We discuss a gravity dual of a charge density wave consisting of a U(l) gauge field and two scalar 
fields in the background of a Schwarzschild-AdS4 black hole together with an antisymmetric field 
(probe limit). Interactions drive the system to a phase transition below a critical temperature. We 
numerically compute the ground states characterized by modulated solutions for the gauge potential 
corresponding to a dynamically generated unidirectional charge density wave in the conformal field 
theory. Signatures of the holographic density waves are retrieved by studying the dynamical response 
to an external electric field. We find that this novel holographic state shares many common features 
with the standard condensed matter version of charge density wave systems. 

PACS numbers: 11.25.Tq, 04.70.Bw, 71.45.Lr, 71.27.+a 



The AdS/CFT correspondence has become a powerful 
tool in studying strongly coupled phenomena in quantum 
field theory using results from a weakly coupled gravity 
background. According to this correspondence principle 
a string theory on asymptotically AdS spacetime is 
dual to a conformal field theory on the boundary. Re- 
cently, following a more phenomenological approach, a 
plethora of condensed matter phenomena have been stu- 
died within this holographic principle, like the properties 
of superfluids and superconductors both conventional 
and unconventional |3[ , Fermi liquid behaviour [4, non- 
linear hydronamics |5[ , quantum phase transitions [fj and 
transport Q- Thus, classical gravity theories have been 
transformed into a laboratory for exploring condensed 
matter physics from a different perspective [8[- 

The most elucidated example of the application of the 
AdS/CFT correspondence to condensed matter physics 
is the holographic superconductor. The gravity dual of 
a superconductor consists of a system with a black hole 
and a charged scalar field coupled to a Maxwell field. The 
coupling induces a negative effective mass for the scalar 
field leading to an instability 0] , which allows to the black 
hole to have scalar hair at temperatures smaller than a 
critical temperature 0. In the boundary dual theory 
this corresponds to a phase transition at a finite chemical 
potential, where the operator dual to the charged scalar 
field condenses. Considering fluctuations of the vector 
potential and calculating the dynamical conductivity, are 
crucial for unveiling the properties of the condensed state. 

Apart from superconductivity, understanding real con- 
densed matter systems demands the introduction of ad- 
ditional ordered states and their holographic interpre- 
tation. The most important of them are undoubt edly 
charge and spin density waves (CDW and SDW) [10|. 
The development of these states corresponds to the spon- 
taneous modulation of the electronic charge and spin den- 
sity, below a critical temperature (T c ). Density waves 
are widely spread in different classes of materials and one 
may distinguish among them, either orbitally pill or Zee- 
man driven pj| field-induced CDWs, confined |13| and 



even unconventional density waves [14| . The latter play a 
special role in condensed matter physics. Characteristic 
representatives are the d-density wave fDDW ) flal and 
most recently its chiral (d+id DDW) version [T(j, both 
of which have been proposed as prominent candidates for 
explaining the still unidentified pseudogap phase of the 
hole-doped cuprates pjj [l6[ and certain 'hidden' orders 
in heavy fermion compounds fl6l. Il7|. Moreover, density 
waves are also significant due to their stron g te ndency to 
compete or coexist with superconductivity jl8l Il9j . 

In this Letter, we put forward a gravity dual of a con- 
ventional charge density wave. Our theory consists of the 
modulus and phase of a complex scalar field, a Maxwell 
field and an antisymmetric field. We numerically ob- 
tain ground states characterized by modulated solutions 
for the scalar potential corresponding to a CDW in the 
boundary dual theory. The solution for the scalar hair 
corresponds to a condensed dual operator which follows 
the usual Bardeen-Cooper-Schrieffer (BCS) dependence 
in the vicinity of T c . Signatures of the holographic CDWs 
are retrieved by studying the collective modes and the 
dynamical response to an external electric perturbation. 
We find, in the dual boundary theory, a gapped phason- 
polariton collective mode p3, HtJ with no Frolich super- 
current [2l| , which clearly shows that our CDW is intrin- 



sically commensurate and non-sliding, a novel property 
stemming from the gravity-imposed pinning of the scalar 
potential to the phason. 

Our starting point for constructing the holographic 
CDW, is the following Lagrangian density in 3+1 dimen- 
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where is the Maxwell gauge field of strength F = dA, 
B^ v is an antisymmetric field of strength H = dB and 
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are auxiliary Stiickelberg fields. The last term equations, 



is a topological interaction term (independent of the me- 
tric). The presence of the Stiickelberg fields is needed 
for the above Lagrangian to be gauge invariant as it can 
be easily checked if the following gauge transformations 
A -> A + d X (1) , -> u;W + B -> B + d X (2) and 
are applied to Eq. ([T]). We shall fix 
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the gauge by choosing = 0, w' 2 ) = 0. Apart from 
the above gauge symmetries, our model is characterized 
by an additional global U(l) symmetry, i9 — > ?9 + a, that 
corresponds to the translational symmetry. This global 
t/(l) symmetry will be spontaneously broken for T < T c 
in the bulk. The concomitant restoration of this global 
symmetry will give rise to the related Nambu-Goldstonc 
mode, the phason (To| as it is called in condensed matter 
physics. One may alternatively understand this global 
symmetry, by unifying the fields $ and i?, into a com- 
plex field as ^ = <I>e lAl1? , where Ai corresponds to the 
charge of the U(l) translational invariance. In a usual 
condensed matter CDW, when translational symmetry is 
completely broken, the phason is massless and the CDW 
is called incommensurate or sliding. The sliding origi- 
nates from the freely propagating phason that gives rise 
to the Frolich supercurrent. In spite of the dissipation- 
less electric charge conduction, U(l) gauge invariance is 
intact and no Meissner effect arises. On the other hand, 
if the phason is gapped then there is a remnant discrete 
translation symmetry that prevents sliding and supresses 
the Frolich conduction. In this case, the CDW is termed 
commensurate or pinned. 

The above model is strongly motivated by the usual 
condensed matter picture, where a CDW corresponds to 
the creation of electron- hole pairs ■ As a consequence, 
the complex scalar order parameter field that describes 
the transition carries zero electric charge in stark contrast 
to the superconducting case where it is charged. There- 
fore, the U(l) gauge invariance of the electric charge will 
be preserved even in the condensed CDW state. More- 
over, the operator that condenses preserves particle num- 
ber. Therefore, we shall be working with a canonical 
ensemble minimizing the free energy F(T,V, N) rather 
than a grand canonical ensemble as is the case with a su- 
perconductor which involves a particle number changing 
order parameter. The chemical potential will be dynam- 
ically generated (unlike with superconductors where one 
obtains dynamics for a fixed chemical potential). 

Despite the phenomenological character of our theory, 
that follows closely the spirit of Rcf . [ 2} , it is well known 
that there exist consistent truncations of M2 string the- 
ory resulting in a supergravity sector with scalar and 
pseudoscalar fields which are neutral under U(l) gauge 
fields (22)]. Such a truncation, although not consistent, 
was used in [23[ to calculate, using the AdS/CFT corre- 
spondence, the momentum relaxation time in a model 
where the translational invariance was spontaneously 
broken due to impurities. 

By varying the metric g^ u , we obtain the Einstein 
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By varying the other fields, A^, B^ v , $ and we obtain 
the field equations, respectively, 
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To set the boundary conditions for the various fields, 
write the metric near the boundary of the asymptotically 
AdS space in terms of Poincare coordinates as 



ds 2 



dr 2 

— + r 2 (-dt 2 + dx 2 + dy 2 ) 



(5) 



where we set L = 1. 

Suppose that A t = V and B xy = -B yx B ry = -B yr 
are the only non- vanishing components of A^ and B^ u , 
respectively We shall also concentrate on static solu- 
tions. In the asymptotic region ([5]), the field equations 
(j3]) become 
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where we have used the notation d r ='. 
are 



The solutions 
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B rv = B (f>) + ( r - 2 ) 

where A± = | ± <J | + m 2 . The constant term in the 
expansion is the source of the corresponding field. No- 
tice that if the source of the B^ field (B^) is spatially 
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dependent, there is an operator in the CFT dual to B xy 
that condenses. Moreover, in B xy , the coefficient of r~ 4 
is arbitrary and represents the vacuum expectation value 
of a dual operator in the CFT. The latter is similar to 
the metric giving rise to the stress energy tensor of the 
boundary CFT. In what follows, we shall concentrate on 
the simple case in which both d x B^ and B^ vanish. 

We wish to solve the field equations in the probe limit 
which is obtained by letting Newton's constant G — > 0. 
In this limit, the Einstein equations © and the equa- 
tions for the B^ v field ([3]) decouple. The latter become 
V p H p ^ 1 ' = 0. These equations can be solved by fixing 
any two-dimensional plane, therefore, by choosing the 
solution 



B ry = _ B % 
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with all other components vanishing, we shall fix the 
(r, y) plane. This also allows us to set B xy = con- 
sistently in this limit on account of ([7]). We stress here 
that the decoupling of the B pv equation in the probe 
limit accommodates a constant source for the field flSJ) 
allowing the development of a (l+l)-dimensional CDW 
on the boundary as we shall see in the following. 

With the above choice of B pu , there is no backreaction 
to the metric and the Einstein equations ([2]) admit the 
planar Schwarzschild AdS black hole 



ds 2 



-f(r)dt 2 + ^-+r 2 (dx 2 + dy 2 ) , /(r)=r 2 -^ 
fir) r 

(9) 

as a solution (L = 1). The metric © asymptotes to the 
AdS metric (|5|) as r —> oo, as expected. 

Having obtained the metric and the B^ v field, the 
other field equations come from the Lagrangian 
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The resulting field equations are independent of G and 
therefore well-defined in the probe limit G — > 0. In fact, 
in this limit, the last term is of the chiral anomaly type 
in a t — x spacctimc, with the fields depending paramet- 
rically on the radial coordinate. This Lagrangian also 
appears in the low energy effective action of the usual 
(1+1)-D CDWs providing further support to the 

phenomenological origin of our theory. 

The direction of the modulation is defined by the non 
zero components of B pv . Therefore, the Lagrangian ap- 
pearing in Eq. (fTU|) shows that the direction of the mod- 
ulation is parallel to the X— axis and it can be used to 
generate any possible profile of modulated charge density 
in the conformal field theory. Our aim is to show that 
the scalar potential A t = V(r, x) asymptotically is of the 
form V(r, x) ~ /i(x) — ^1 + . . with /z(x) ~ cos{kx) and 
p(x) ~ cos(fca;). According to the AdS/CFT dictionary 
these two terms of this expansion will correspond in the 



dual boundary theory to a chemical potential fJ-(x) and a 
charge density p(x), which are modulated. Consequently, 
all the fields will be x— dependent. The chemical poten- 
tial will be dynamically generated. Thus, we shall be 
working with a canonical ensemble of fixed particle num- 
ber. This is possible in our case because the operator 
that condenses docs not change the particle number. 

Varying the Lagrangian (|10[) we obtain the field equa- 
tions 
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which are also obtained from (|3]) with the background 
choices ((8]) and (0. In the first field equation there is a 
direct coupling of the scalar potential V with the phase 
This has important consequences as we will discuss in 
the following. 

Suppose that the .T-direction has length L x . Then 
assuming periodic boundary conditions, the minimum 
wavevector is 
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Taking Fourier transforms, the form of the field equations 
suggests that it is consistent to truncate the fields by 
including (2n + l)/c-modes for •d and V and 2nk modes 
for $, where n € Z. Thus 
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We shall consider only the contributions of the leading 
(n = 0) terms #±fc, V±k and $o- Higher overtones can 
be included successively. At each step, we need to add 
equations for the new overtones which also introduce cor- 
rections to the lower Fourier modes already calculated in 
the previous step. Evidently, at each step the complexity 
of the system of coupled equations for the overtones of 
the three fields increases. We plan on carrying out this 
systematic procedure elsewhere. Ignoring overtones, the 
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At the boundary, 
where A = 1 , 2 



We shall solve the above equations in the case where the 
mass of the scalar field is m 2 = —2, which is above the 
Breitenlohner-Freedman bound [25[ . We shall also fix the 
parameters Ai =A2 = 1. In order to find the boundary 
conditions at the horizon we take the near horizon limit 
r — > rhi where f(r) ~ — Hi), and extracting the 

parts that become singular in the above equations we find 
that V± k (r h ) = #± k (r h ) = 0, %(r h ) = m 2 $ (r h )/3r h 
and d' ±k {r h ) = ^ikV^r H ) / $r h f 

we have asymptotically <&o(r) 

in general, on account of ([7J. In this work we focus on 
A = 1 leaving the other case to a future study. 

The fields ■& and V are also normalizable and behave 
asymptotically as in ([7} with A = 1. The coefficients 
in their respective expansions have a sinusoidal depen- 
dence on x. Thus we obtain in the dual boundary the- 
ory a single-mode CDW with a dynamically generated 
charge density of the form p{x) = V^ 1 ' ~ cos(fcx). This 
case is important because a lot of materials exhibit cither 
a single-mode CDW or only one dominant wave-vector. 
Moreover, the resulting model is numerically tractable 
and can provide direct insight on the phase transition 
and the response of this novel holographic state. 

Having determined the behaviour of the fields both 
at the horizon and the boundary, we solve the five field 
equations ([T4^) - ([T6l numerically setting k = 1. We clearly 
obtain a second-order phase transition with a modulated 
scalar potential. Notice that a spatially modulated phase 
transition has been also observed in a five-dimensional 
Maxwell theory with a Chern-Simons term (2(|. In Figfl] 
we plot the temperature dependence of the condensate 
(0i ) near the critical temperature. By fitting with |T — 
Tel* 3 , we find that the critical exponent for the transition 
is 0.5, i.e. of the BCS type. 

The critical temperature is controlled by the value of 
the wave-vector k. Its dependence on k can be deduced 
from the scaling symmetry of the system under the trans- 
formation t — > (t, x — > (x, y — »• Cy, r — > r/£, V — > V/(. 
It follows that k/rh is scale- invariant, therefore 



T c cx k 



(17) 



We have checked the validity of this conclusion numeri- 
cally. Eq. (fTTf should be contrasted with the case of a 
holographic superconductor in which T c cx yfp. 




FIG. 1: Temperature dependence of the condensate (solid 
line). The dashed line is the BCS fit to the numerical values 
near T c . We find (Oi) w 8.5T C (1 - T/T c ) 1/2 near T -> T c . 



Also, notice that when the condensate (0i) is zero, 
i.e. for temperatures above the critical temperature, 
the modulated chemical potential and the charge density 
vanish, and that they become non zero as soon as the 
condensate becomes non zero, i.e. when the temperature 
is lowered below T c . Therefore, the modulated chemical 
potential and the charge density are spontaneously gene- 
rated and do not constitute fixed parameters of control- 
ling T c , contrary to what happens in holographic super- 
conductors 0] with their homogeneous analogs fi, p. We 
also have to remark that the numerical solution of these 
five equations becomes unstable at low temperature. In 
fact, the fields diverge for T/T c < 0.4. This numerical 
instability indicates that the probe limit breaks down at 
low temperatures and backreaction effects on the bulk 
metric become significant. A further analysis away from 
the probe limit is certainly needed. 

We proceed with examining the collective excitations 
and response in the dual boundary theory by considering 
fluctuations of the electromagnetic and phason fields of 
the form A x (r,x)^> A^^e 1 ^-^, V(r,x)->V(r,x) + 

V Ut q(r)e i( - qx - ut ) and ^(r,x)^(r,x) + d^^e^-^. 
In the Lorentz gauge, we obtain the equations 
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= 0. 



(19) 
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(20) 



In order to get an insight of the response of the our system 
to these fluctuations we take the limit 0. After 
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setting z = 1/r and Fourier transforming z — > p, we 
obtain 

(lu 2 -q 2 - p 2 ) ? Ui ,(p) + igK,,,(p) + = ,(21) 

+*g(Oi) 2 ^, 9 (p) + {lo 2 -q 2 - p 2 ) V u<g (p) = ,(22) 
-iu;(0i) 2 ^, 9 (p) + - <1 2 - V 2 ) A u . q (p) = .(23) 

The above homogeneous system of equations defines com- 
pletely the dynamics of the fluctuations through simply 
setting its determinant equal to zero. Hence, we get 

lu 2 = q 2 + p 2 , u? = q 2 + p 2 + lu 2 ± iu g y^Tv , (24) 

where we have defined ^/2uj g = (Oi). Each Fourier 
p— mode of "&u>,q, Vu,q and A UtQ will give a solution of the 
form e l r. with p = p(cj,q). Near infinity e 1 ^ ~ 1 + if, 
and according to the AdS/CFT correspondence, the first 
term with p = corresponds to a source and the second 
to a current in the boundary theory. This implies that 
the dispersions of the collective modes must be calculated 
for p = 0, while for p ^ wc obtain the corresponding 
current that controls the response of the system such as 
the conductivity. 

For p = 0, we obtain from Eq. ((24)) three energy dis- 
persions uj = q, lu = q, lu = \J {Oi) 2 + q 2 . The first two, 
correspond to masslcss photonic-like modes that we an- 
ticipated to find since gauge invariance still persists. On 
the contrary, the last mode has a mass for q = equal 
to (Oi ) and basically corresponds to a gapped phason- 
like mode which originates from phason-gauge coupling 
[20l | . The emergence of the gap hints to the pinning of the 
CDW [Il|,[24|. As a matter of fact, this pinning has an in- 
trinsic origin, because when the phase transition occurs, 
both the scalar potential and the phason field become fi- 
nite and modulated by the same wave- vectors ±fc. Their 
existence is inseparable and in stark contrast to the holo- 
graphic superconductor 0, the scalar potential is zero 
above T c . The interaction term V (x)d x 'd(x) in Eq. (flQl) . 
demands that if V(x) ~ cos(fca;) then $(x) ~ sin(foc). 
Their relative phase is 'locked' at \ and the CDW is 
commensurate . 

The dynamical conductivity (<j = 0) is dominated by 
strongly hybridized J? w ,o(jP)j A u fi(jp) fluctuations and is 
defined from Ohm's law, according to the AdS/CFT 

correspondence, as cr(w) = A^' /iwA]y , where A^ 
and A^q are determined from the asymptotic expansion 

~ ' ~(o) 

A U fi(r) — A u o H ^ + • • • . The Fourier expansion 

Auj.o{r) = J2 P Aj,o(p)e 1 ^, implies that asymptotically 

AE/o = iY.pPAuflip) and A^ Q = £ p A Wj0 (p). Fur- 
thermore, we obtain from Eq. (|24l) . for q = in the 
second equation, that p±.± = ±^/lu(lu ± (Oi)), which 
correspond to two different branches for propagation di- 
rection (first sign). By choosing only the ingoing con- 
tributions p—,± we obtain the dynamical conductivi- 
ty <t(u) = T,e=± £ ^P 1 where c e = Aufi(P-,e)/A { °l. 
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FIG. 2: The numerically calculated real part of the conduc- 
tivity _Re[<r(a>)], versus normalized frequency with the con- 
densate u/(Oi) (left) and temperature oj/T (right). In both 
plots, we clearly observe a 'dip' that arises from the CDW for- 
mation and softens with T — > T c since (Oi) — > 0. At T = T c 
we retrieve the normal state conductivity Re[a(uj)] — 1. 

The factors c e will be determined by demanding that 
lim w ^oo [f(^) — 1] — > faster than 1/w, in order for the 
Kramers-Kronig relations to hold. The condition implies 
c e = i which we have verified numerically that also holds 
for finite temperatures. Kramers-Kronig relations require 
the fulfillment of the Ferrell-Glover-Tinkham (FGT) sum 
rule, dictating that / dw Re[cr(cj)] is a constant indepen- 
dent of temperature. In the superconducting case 0, the 
FGT sum rule demands the presence of S(lu) in Rc[ct(o;)], 
giving rise to a supcrcurrent. In our case, this rule is 
satisfied exactly, without the need of S(lu). The latter re- 
flects the absence of the Frolich supercurrent, which may 
be attributed to the commensurate nature of the CDW. 

We now present numerical results for the dynamical 
conductivity in the vicinity of T c . In this limit, we are 
left to solve the system of the coupled equations Eq. 
(fl~8)) and (|20|) . The stable ingoing boundary conditions 

at the horizon read A U)0 (r h ) oc £± f-W^ ±r ^/ 3r » 

and d u , (r h ) oc £± i(±r h $)- 1 f-hfitp^^*)/*-* . On 
the first panel of Fig[2] we plot the real conductivity ver- 
sus u/(Oi) for different temperatures, while on the sec- 
ond, the same data are shown versus lo/T. Notice how 
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the phason-polariton interplay results in a 'dip' structure 
which gets sharper as temperature decreases. Moreover, 
this feature moves to higher frequencies as the conden- 
sate value increases for lower temperatures. Upon raising 
the temperature, the 'dip' CDW fingerprint washes out, 
disappearing totally at T = T c where the normal state is 
restored. 

In conclusion, we have proposed a holographic charge 
density wave. The dual boundary theory has the charac- 
teristic features of a commensurate CDW, with a gapped 
phason mode and the absence of a Frolich supercurrent. 
The commensurate behaviour due to the locking of the 
relative phase of the scalar potential and the phase field is 
an intrinsic property of the system and may be attributed 
to the presence of the black hole in the bulk. Moreover, 
the dynamical conductivity shows phason-polariton mi- 
xing. The properties of the holographic charge density 
wave certainly need further investigation, as for example 



the possibility of a multi-fc CDW, the effect of impuri- 
ties, the response to an external magnetic field and its 
competition to superconductivity. Moreover, it would be 
interesting to examine this model away from the probe 
limit and whether this phenomenological model can be 
obtained from a consistent truncation of string theory. 
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